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Abstract 


The occurrence of land subsidence phenomenon and its potential hazards in the plains of Iran 
due to the water crisis and drought period has grown significantly in recent years. In this study, 
the zoning of land subsidence hazard in the Hashtgerd plain was discussed and 19 criteria were 
selected as factors that influence the subsidence. The mentioned layers were prepared in GIS, 
weighted based on the best-worst method (BWM) and integrated using the Weighted Overlay 
Index (WOI). The results of BWM method showed that groundwater abstraction (0.219), the 
type of geological formation (0.157), the decrease of groundwater level (0.079), and 
groundwater depth (0.078) are important factors on the potential of subsidence hazard in the 
study area. Moreover, in order to evaluate the results of this model, the ROC curve was used, 
which has an accuracy of 90%. The results showed that land subsidence hazard of 10.66% of the 
study area was in very low category, 38.51% in low category, 31.49% in medium category, 
11.66% in high category and 7.69% is in the very high category. According to the results, there 
are areas with a high probability of subsidence in the central part of Hashtgerd Plain, which 
requires continuous evaluation, controling, and monitoring of the criteria that affect the 
situation. 
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جغرافیا و مخاطرات محیطی. سال یازدهم. islo‏ چهل و چهارم؛ زمستان ۱۶0۱ صص ۱۲۷-۱۶۸ 
مقاله پژوهشی 
ویژه‌نامه (چالش جهانی فرونشست زمین: مدیریت بحران یا بحران مدیریت) 


پهنه‌بندی خطر وقوع فرونشست زمین در دشت هشتگرد بر اساس رویکرد تلفیقی تصمیم گیری چند معیار: 
WOI-BWM‏ 
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حکیده 


وقوع پدیده فرونشست زمین و خطرات احتمالی آن در دشت‌های ایران به علت بروز بحران آب و 
خشکسالی» طی سالیان اغیر رشد چشمگیری داشته است. در این مطالعه به پهنه‌بندی خطر وقوع 
فرونشست در دشت هشتگرد پرداحته و نوزده معیار شامل میزان افت آب زیرزمینی برداشت آب 
زیرزمینی. ضخامت آبخوان. ضخامت آبرفت. تغذیه آبخوان, تراکم چاه‌ها؛ تراکم زهکشی. عمق سطح آب 
زیرزمینی» نوع سازند زمین‌شناسی» عمق سنگ‌بستر متوسط بارش سالانه. متوسط دمای سالانه نوع اقلیم 
سرزمین کاربری کشاورزی» کاربری شهری, کاربری صنعتی فاصله از آبراهه» فاصله از جاده و فاصله از 
گسل به‌عنوان عوامل موثر بر وقوع فرونشست منطقه انتخاب شدند. clea‏ مذکور در سامانه اطلاعات 


جغرافیایی تهیه بر اساس مدل THE‏ بدترین Öja (BWM)‏ دهی و با استفاده از شاخحص هم‌پوشانی 
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۱۳/۸ جغرافیا و مخاطرات محیطی شمارة چهارم 


وزین (WIO)‏ ادغام شدند. نتایج روش BWM‏ نشان داد عوامل برداشت از آب زیرزمینی TNA)‏ نوع 
سازند زمین‌شناسی (۰/۱۵۷ افت سطح آب زیرزمینی (۰/۰۷۹) و عمق آب زیرزمینی (۰/۰۷۸) اهمیت 
بیشتری بر پتانسیل خطر فرونشست در منطقه موردمطالعه دارند. همچنین به‌منظور ارزیابی نتایج این مدل 
از منحنی (ROC)‏ استفاده شد که دارای دقتی معادل ٩۰‏ درصد است. نتایج نشان داد که ۱۰/17 ۳۸/۵۱ 
۹ ۱۱/۷ و ۷/۹۹ درصد از مساحت دشت هشتگرد به ترتیب به طبقه خطر خیلی کې کم متوسط 
پرخطر و خیلی پرخطر اختصاص atl,‏ است که بر این اساس مناطقی با احتمال فرونشست زياد در نیمه 
مرکزی دشت هشتگرد قرار دارند که مستلزم ارزیابی نظارت و GUL‏ مستمر عوامل مؤثر در وقوع پدیده 
فرونشست توسط مدیران و برنامه‌ریزان شهری است. 
کلیدواژه‌ها: فرونشست زمین, تغذیه آبخوان. دشت هشتگرد. مدل BWM‏ 


۱- مقدمه 


در سال‌های اخیر» یک خطر زمین‌شناسی که در بسیاری از شهرهای جهان گسترش یافته است. مسأله فرونشست 
زمین است (جن ' و همکاران. ۲۰۱۹). فرونشست زمین» همان‌طور که توسط یونسکو تعریف شده است. شامل 
نشست یا نشست تدریجی سطح زمین به سمت پایین است که ممکن است با جابه‌جایی‌های عمودی با بردار اندک 
افقی نیز همراه باشد (یونسکو " ۲۰۱۸). فرونشست زمین یکی از پیامدهای استفاده بی‌رویه و سوء مدیریت از منابع 
آبی است که به دلیل فشردگی خاک در اثر کاهش سطح آب‌های زیرزمینی رخ می‌دهد (پاچکو" و همکاران, ۲۰۰3) 
که باعث آسیب به محیط طبیعی و حتی خسارات اقتصادی می‌شود (هو* و همکاران. ۲۰۰۹). در فعالیت‌های انسانی 
این پدیده متأثر از فعالیت‌هایی مانند برداشت بیش‌ازحد آب از سفره‌های زیرزمینی Sly)‏ و همکاران ۲۰۱۹ 
ندیری" و همکاران» ۲۰۱۸)؛ زهکشی خاک‌های آلی (توسی " و همکاران ۲۰۱۳) استخراج سوخت‌های فسیلی و 
استخراج زیرسطحی همچون معدنکاری (گالوی ۲۰۱7؛ افزایش حفاری چاه‌های بهره‌برداری (ژانگ* و همکاران؛ 
۹ و پدیده‌های طبیعی زمین‌شناسی همچون حرکات تکتونیکی گسل‌ها (سوپاتا " و همکاران. ۲۰۲۰؛ پرادهان ۲ 
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سال یازدهم پهنه بندی خطر وقوع فرونشست زمین در دشت ... ۱۳۹ 


و همکاران» ۲۰۱۶) و کارستی شدن به دلیل انحلال سنگ‌آهک (دزیر" و همکاران: ۲۰۱۸) است. فرونشست اگرچه 
پیامد فشردگی سطوح زیرین خاک است لیکن خود معرف حدوث و پیامد ناگوار زیست‌محیطی در افق آتی است. 
فشردگی سطوح زیرین و افت عمق ایستابی» کاهش حجم و فضای نگهداشت آب را به دنبال خواهد داشت. به بیان 
ساده‌تر کاهش مخزن آب زیرزمینی پیامد منطقی رخداد پدیده فرونشست است (شریفی‌کیا و همکاران ۱۳۹۶) که 
می‌تواند خسارت‌های شدیدی به تأسیسات انسان‌ساخت (ساختمان)» باند فرودگاه‌هاء پل‌هاء تونل‌هاء خیابان‌ها؛ 
خطوط حمل‌ونقل ریلی و جاده‌ای» کاهش حاصلخیزی زمین‌های کشاورزی L)‏ فشرده شدن و از بین رفتن تخلخل 
موجود در (SEE‏ اختلال در الگوی جریان‌ها و تداخل آب‌شور و شیرین و تخریب مستحدثات زیرسطحی (قنات‌ها؛ 
شبکه سوخت و آب‌رسانی» تونل‌های شهری و معادن) را به دنبال داشته باشد (شفیعی و همکاران ۱۳۹۸). 
ازآنجایی که اولین گام در مطالعه این پدیده شناخت گستره فضایی آن و اندازه گیری هر چه دقیق‌تر میزان فرونشست 
است؛ بنابراین با شناخت مشخصات مکانی و زمانی این پدیده. امکان ارائه و توسعه مدل منطقه‌ای این پدیده میسر 
شده (مهرابی و غضنفری‌پور ۱۳۹۸) و از این طریق می‌توان زیرساخت مناسبی از منابع داده و اطلاعاتی را برای 
مدیران و دست‌اندرکاران در راستای تعدیل» JES‏ و پیشگیری از فاجعه فراهم آورد (شریفی‌کیا و همکاران» ۱۳۹۲). 
مطالعات متعددی برای ارزیابی فرونشست با استفاده از سه رویکرد: الف) مدل‌های ریاضی با استفاده از مدل‌های 
جریان آب زیرزمینی و مدل‌های فرونشست برای بررسی اثرات کاهش سطح آب منجر به فرونشست؛ ب) 
تکنیک‌های سنجش‌ازدور برای نظارت بر فرونشست زمین با قابلیت بررسی روابط احتمالی با رویکرد 
هیدروژئولوژیکی؛ ج) استفاده از تکنیک‌های GIS‏ برای نقشه‌برداری مناطق آسیب‌پذیر در برابر فرونشست انجام شده 


است (صادق فام" و همکاران. ۲۰۲۰). در جدول (۱) سابقه‌ای از آن ارائه شده است. 


جدول ۱- مرور منابعی از مطالعات مرتبط با پدیده فرونشست 


طبقه نویسندگان داده‌های مورداستفاده موقعیت نتایج کلیدی 
Sse‏ تغییر شکل زمین. خحصوصیات ۳ s‏ 
مدنی " و همکاران US ob BS‏ همبستگی بین فرونشست زمین, برداشت اب 
۱ خاک زیرین و رژیم اب‌های 1 ۱ ۱ 
زیرزمینی 
ری 2 7 ۲ 
فانت* و همکاران | جریان آب‌های سطحی, نقشه‌های کالیفرنیا ارزیابی استراتژی‌های مدیریتی برای کاهش اثرات 
(Ye)‏ کاربری زمین و داده‌های آب‌وهوا (آمریکا) نامطلوب فرونشست صورت گرفت. 
Desir‏ 1 
Sadeghfam‏ 2 
Modoni‏ 3 


4 Faunt 


۱۳۰ جغرافیا و مخاطرات محیطی شماره چهارم 
aib‏ نویسندگان داده‌های مورداستفاده موقعیت نتایج کلیدی 
خحصوصیات فیزیکی, اقلیم» تغذیه دره 
شرستها" و , پهنه‌بندی خطر فرونشست. مکان و سطوح بالقوه خطر 
اب‌های زیرزمینی تراکم‌پذیری و | کاتماندو 2 
همکاران (۲۰۱۷) ۱ را برای مدیریت پیشگیری از خطر برجسته می‌کند. 
نوسان سطح آب‌های زیرزمینی (Sk)‏ 
7 میزان فرونشست دارای تغییرات فضا و زماد ابیت 
آبیدین" و همکاران تصاویر راداری INSAR‏ جاکارتا ات 
که ناشی از استفاده بیش‌ازحد از آب‌های زیرزمینی» با 
سی از بیس‌ازحد از اآب‌های زیرزمینی بار 
)۰10( و GPS‏ (اندونزی) ~ 
سازه‌ها و فعالیت‌های زمین‌شناسی در منطقه است. 
منطقه به دلیل نداشتن موانع ژئومورفولوژی» در طی 
سال‌های اخیر با روند توسعه ساخت‌وسازها؛ توسعه 
ابراهیمی و پاکدشت 
تصاویر ستتینل ۱ فعالیت‌های صنعتی و کشاورزی مواجه شده و همین 
همکاران (۱۳۹۹) (ایران) = 
امر سبب افزایش بهره‌برداری و پدیده فرونشست شده 
است. 
مدل‌ها L‏ 
led‏ بر ah‏ ارتفاع شيب» جهت؛ 
Ga‏ سنگ‌شناسی, فاصله از گسل 
پرادهان "و دره کینتا ‏ | نتایج نشان داد که مدل ۲8 در پیش‌بینی نسبت به مدل 
: فاصله از رودخانه. Sle NDVI‏ ۱ 
همکاران (۲۰۱۶) (مالزی) تابع باور شواهد دقت کمتری دارد. 
قدرت = Ob‏ رودخانه» رطوبت. 
توپوگرافی 
فاصله از قنات. کاربری اراضی؛ ۱ ۱ 
bol ۱ ۳ ۱ :‏ بیان کردند برداشت GR!‏ زیرزمینی و به دنبال 
رحمتی و فاصله از زمین‌های جنگلی. کاشمر 5 E‏ 
۳ ۳ ان سنگ‌شناسی و فاصله از رودخانه‌ها مهم‌ترین عوامل 
همکاران (۲۰۱۹) | سنگ‌شناسی فاصله از گسل‌هاء (ایران) 


در فرونشست هستند. 


سطح آب زیرزمینی 


در حال حاضر این پدیده به همراه پیامدهای زیست‌محیطی آن در اکثر دشت‌های دنا به چشم می‌خورد (رهنما 


و میرائی" TV‏ در ایران نیز فرونشست سابقه طولانی دارد. در واقع با توجه به وضعیت نیمه‌عشک بسیاری از 


مناطق ایران و نیاز به منابع cal‏ میزان بهره‌برداری از آب‌های زیرزمینی بسیار بیشتر از ميزان تغذیه آن‌ها است» به 


(ابراهیمی و همکاران. ۳۹۹۵ در این میان» یکی از مناطقی که با مخاطره فرونشست مواجه است» دشت هشتگرد در 


غرب استان تهران است که با توجه به اهمیت موضوع در این تحقیق به ارزیابی میزان فرونشست آن پرداخته شده 


است. بررسی‌های سازمان زمین‌شناسی و اکتشافات معدنی Oly!‏ (۱۳۸۹) با استفاده از روش تداخحل سنجی ماهواره‌ای 


1 Shrestha 

2 Abidin 

3 Pradhan 

4 Rahmati 

5 Rahnema & Mirassi 


سال یازدهم پهنه بندی خطر وقوع فرونشست زمین در دشت ... ۱۳۱ 


امواج رادار نشان داده است که در پهنه‌ای گسترده از این دشت. فرونشست زمین با نرخ بیشینه ۱٩‏ و نرخ میانگین 
سطح آب زیرزمینی مواجه بوده است. همچنین گمان می‌رود که در آینده نیز با توجه به تغییرات اقلیمی و با فرض 
ادامه یافتن شرایط موجود در برداشت از منابع آب زیرزمینی» کاهش سطح آب‌های زیرزمینی و بالطبع آن پدیده 
فرونشست با شدت و گستردگی بیشتری ادامه یابد. لذا تهیه نقشه پهنه‌بندی خطر فرونشست برای حفاظت از منابع 
تعیین مهم‌ترین عوامل موثر در وقوع فرونشست با استفاده از روش وزن‌دهی BWM‏ و ارزیابی صحت این مدل در 


دشت هشتگرد انیت 


۲- مواد و روش‌ها 

۱-۲- موقعیت منطقه موردمطالعه 

دشت هشتگرد از دیدگاه تقسیمات کشوری در استان البرز قرار دارد. مساحت این حوضه در حدود ۱۱2۸/۸ 
کیلومترمربع است که در عرض جغرافیایی "10 EV!‏ ۳۵۹ تا ۰۳۰۵۳ ۳۹ شمالی و طول جغرافیایی ۰۵۳ ۲۹۲ *۵۰ تا 
۳ 087 ۵۰3 شرقی قرار دارد. آب‌وهوای منطقه در بخش‌های شمالی» نیمه مرطوب است و به‌تدریج به‌سوی جنوب 
دشت با کاهش ارتفاع به نیمه‌حشک تبدیل می‌شود (شکل ۱). بررسی‌های زمین‌شناسی منطقه برای شناخت پدیده 
فرونشست ضروری است. به‌طورکلی از دیدگاه هدایت هیدرولیکی. سازندهای زمین‌شناسی این محدوده به ٤‏ گروه 
اصلی تقسیم می‌شوند: الف) سازندهایی با نفوذپذیری بسیار کم؛ ب) سازندهایی با نفوذپذیری کم؛ پ) سازندهایی 
با نفوذپذیری متوسط و ت) سازندهایی با نفوذپذیری زیاد تا خیلی زیاد. همچنین ساختار زمین‌ساختی ناحبه 
موردمطالعه از فعالیت‌های زمین‌ساختی البرز تأثیر پذیرفته است. به همین دلیل» بیشتر چین‌خوردگی‌ها و گسل‌ها؛ 
روند خاوری- باختری دارند. از گسل‌های مهم این منطقه می‌توان گسل آبیک. گسل جنوب طالقان و گسل ولیان را 


نام برد (سازمان زمین‌شناسی و اکتشافات معدنی YAY cil pl‏ 
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۲-۲- روش انجام پژوهش 

پژوهش حاضر در پنج مرحله کلی به انجام رسیده است که شامل: ۱) شناسایی متغیرهای محیطی تأثیرگذان ۲) 
آماده‌سازی لایه‌های اطلاعاتی مربوط به متغیرهای محبطی موژثر بر وقوع فرونشست زمین؛ ۳) بررسی اهمیت 
متغیرهای محیطی بر اساس روش بهترین- بدترین E (BWM)‏ اجرای مدل شاخص همپوشانی وزین و تهیه نقشه 
پهنه‌بندی فرونشست زمین و ۵) اعتبارسنجی مدل با استفاده از منحنی تشخیص عملکرد نسبی (ROC)‏ است. مراحل 
کلی این پژوهش در شکل نمایش داده شده است (شکل AY‏ 
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وزن‌دهی براساس 
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شکل ۲- مراحل انجام فرآیند پهنه‌بندی فرونشست در دشت هشتگرد 


۱۳ جغرافیا و مخاطرات محیطی شمارة چهارم 


۱-۲-۲- تعیین معیارهای مؤثر بر فرونشست 

در این پژوهش, با توجه به شرایط دشت هشتگرد yy‏ منابع و همچنین نظر افراد خبره معیارهای تأثیرگذار در 
پهنه‌بندی فرونشست زمین در ۲ بعد فیزیکی- شیمیایی و اقتصادی- اجتماعیء ۵ گروه معیارهای هیدرولوژی؛ 
زمین‌شناسیء هوا و اقلیم. کاربری اراضی و حرایم و ۱۹ زیرمعیار شامل میزان افت آب زیرزمینی؛ برداشت آب 
زیرزمینی. ضخامت آبخوان» ضخامت آبرفت. تغذیه آبخوان, تراکم چاه‌ها؛ تراکم زهکشی. عمق سطح آب زیرزمینی؛ 
نوع سازند زمین‌شناسی. عمق سنگ‌بستر میزان بارش سالانه متوسط دمای سالانه نوع اقلیم سرزمین» SAMS‏ 
کشاورزی» کاربری شهری کاربری صنعتی. فاصله از آبراهه‌هاء فاصله از جاده‌ها و فاصله از گسل‌ها تدوین شدند. 
عابدینی و همکاران (VEN)‏ در پژوهشی بیان کردند که افت سطح آب. فاصله از رودخانه و لیتولوژی از مهم‌ترین 
عوامل موْثر در وقوع فرونشست زمین در دشت اردبیل است. محمدی " و همکاران (۲۰۱۹) برداشت از آب‌های 
زیرزمینی» ساختار زمین‌شناسی. فاصله از گسل و درصد شیب را از عوامل موثر در فرونشست بیان کردند. 

۲۷-۲۷-۲- آماده‌سازی معیارها 

با تعیین مجموعه‌ای از معیارها؛ نیاز است که هر معیار به‌صورت یک لايه نقشه در پایگاه داده‌های مبتنی بر GIS‏ 
نشان داده شوند. در ذیل به تشریح معیارها خواهیم پرداخت: 

- فاصله از آبراهه‌ها: آبراهه‌ها با افزایش فشار آب منفذی یکی از عوامل مهم در القاء ناپایداری در منطقه است 

(رنجگر" و همکاران» ۲۰۲۱؛ قربان‌زاده" و همکاران, ۲۰۲۰). لذا نقشه آبراهه‌ها از نقشه توپوگرافی ۱:۵۰۰۰۰ 
منطقه موردمطالعه استخراج گردید و سپس با استفاده از ابزار Distance‏ در نرم‌افزار ArcGIS‏ نقشه فاصله از آبراهه 
تهیه گردید که مقادیر پیوسته آن بین صفر تا 1۸۲۳ متر است. 

فاصله از گسل: گسل‌ها نقش بسیار مهمی در شکستگی و خوردشدگی واحدهای زمین‌شناسی و حرکت‌های 
توده‌ای دارند (نگهبان و همکاران» ۱۶۰۰). لرزش ناشی از ایجاد گسل می‌تواند باعث تشکیل یک تراکم غیر الاستیک 
به دلیل افزایش تنش موثر در SLE‏ و ایجاد فرونشست شود (گالووی" و همکاران» ۹۹۹ ناوارو هرناندز" و 
همکاران؛ ۲۰۲۰). در این پژوهش. گسل‌های منطقه از نقشه زمین‌شناسی با مقیاس ۱:۱۰۰۰۰۰ استخراج گردید و 
سپس با استفاده از ابزار Distance‏ در نرم‌افزار ArcGIS‏ نقشه فاصله از گسل تهیه گردید. 

میزان برداشت آب و افت سطح آب‌های زیرزمینی: یکی از دلایل فرونشست زمین مربوط به عکس‌العمل طبیعی 
آن در مقابل بهره‌برداری بیش ازاندازه و بالطبع آن افت سطح آب‌های زیرزمینی به دلیل کاهش فشار هیدرولیکی و 
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افزایش تنش بر رسوبات است (هانگ ' و همکاران» ۲۰۱۲؛ هو" و همکاران. ۲۰۰۹). در این پژوهش به‌منظور تهیه 
نقشه میزان افت سطح آب‌های زیرزمینی محدوده مطالعاتی از اطلاعات تغیبرات دوره ۱۵ ساله عمق آب از شرکت 
آب منطقه‌ای استان البرز (۱۳۹۹) استفاده شده است و سپس با استفاده از روش درون‌یابی IDW‏ نقشه پیوسته ميزان 
افت سطح آب زیرزمینی تهیه شد. همچنین» میزان برداشت آب زیرزمینی با استفاده از داده‌های تخلیه چاه‌های 
بهره‌برداری تهیه گردید و سپس با استفاده از روش درونیابی IDW‏ نقشه پیوسته برداشت آب زیرزمینی تهیه شد. در 
نقشه پیوسته مقادیر آن بین ۷ تا ۱۷۵۲ مترمکعب است. 

تغذیه: تغذیه یکی از عوامل هیدرولوژیکی مؤثر در فرونشست است. هرچه میزان تغذیه بیشتر باشد» فشار 
هیدرولیکی بیشتر شده و فاصله بین دانه‌ها افزایش می‌یابد که در نتیجه تنش مؤثر کاهش یافته و احتمال وقوع 
فرونشست کمتر می‌شود (صادق‌فام "9 همکاران ۲۰۲۰؛ ندیری* و همکاران. ۲۰۱۸). در این پژوهش. جهت تهیه لایه 
تغذیه از روش پیسکوپو (پیسکوپو ۲۰۰۱) استفاده شد. در روش پیسکوپو لایه شیب زمین» بارندگی و نفوذپذیری 
سطحی خاک (عمق ۰/۵ تا ۲ متری SE‏ سطحی) به a‏ رستری تبدیل و با وزن‌دهی و رويهم‌گذاری لایه‌هاء نقشه 
پیوسته تغذیه تهیه شد که مقادیر آن بین ۰/۳۹ تا یک است. 

ضخامت آبرفت: این ضخامت شامل فاصله بین سطح زمین تا سنگ کف است. در حالت کلی مشاهده می‌شود 
که نشست بیشتر در plea‏ اتفاق می‌افتد که ضخامت آبرفت بیشتر باشد زیرا اگر ضخامت آبرفت منطقه بیشتر شود 
آب زیادی می‌توان استخراج کرد. برای ad‏ این لایه از داده‌های ژئوفیزیک سطحی و چاهپیمایی که از شرکت آب 
منطقه‌ای استان البرز (۱۳۹۹) SEN‏ شده بود میزان ضخامت آبرفت در هر چاه منطقه محاسبه گردید و سپس این نقاط 
بر اساس روش DW‏ درون‌یابی شد و به فایل رستری تبدیل گردید و بدین ترتیب لایه پیوسته ضخامت آبخوان 
به‌دست آمد که مقادیر آن بین ۰" تا ۳۲۹ متر است. 

پارامترهای اقلیمی: پارامترهای اقلیمی به خصوص دما و بارش به دلیل تأثیری که بر آب‌های زیرزمینی دارد؛ یکی 
از عوامل مؤثر در پدیده فرونشست است (ژنگ" و همکاران. ۲۰۱۷؛ تفرشی" و همکاران ۲۰۱۹). لذا در این 
پژوهش. داده‌های جدولی همانند موقعیت و اطلاعات ایستگاه‌های هواشناسی مربوط به میانگین بارندگی سالیانه و 
میانگین cles‏ سالیانه منطقه موردمطالعه به‌صورت نقطه به محیط نرم‌افزار ArcGIS‏ وارد گردیده و سپس با استفاده از 


روش‌های درون‌یابی IDW‏ به نقشه هم‌باران سالبانه و هم‌دمای alle‏ تبدیل گردید. مقادیر نقشه oble‏ بین ۱۹۹ تا 
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۰ میلی‌متر در سال و مقادیر هم‌دما بین صفر تا ۱۶ درجه سانتی‌گراد است. همچنین» فرونشست معمولاً در 
اقلیم‌های خشک که هم برداشت از آب زیرزمینی بیشتر است و هم تغذیه سفره آب زیرزمینی ضعیف‌تر است اتفاق 
می‌افند. نقشه نوع adil‏ در منطقه موردمطالعه از سازمان هواشناسی استان البرز دریافت گردید که بر اساس آن اقلیم 
منطقه دارای " کلاس است که کلاس اول معرف خشک‌ترین و کلاس ششم معرف مرطوب‌ترین کلاس اقلیمی است. 

سازندهای زمین‌شناسی: سنگ‌ها به‌واسطه تفاوت در جنس رسوبات تشکیل‌دهنده و شرایط و دوران شکل‌گیری» 
مقاومت‌های مختلفی را در مقابل نیروهای خارجی از خود نشان می‌دهند (کانفورتی ! و همکاران NONE‏ مختاری و 
عابدیان ؛ ۲۰۱۹). به همین منظور, نقشه رقومی سنگ‌شناسی منطقه با مقیاس ۱:۱۰۰۰۰۰ از سازمان زمین‌شناسی و 
اکتشافات معدنی (۱۳۹۹) دریافت گردید و خصوصیات واحدهای سنگ‌شناسی استخراج گردید (جدول (Y‏ منطقه 


مج 


هشت‌گرد دارای ۱۹ کلاس سازند زمین‌شناسی ات که سازندها بر اساس جنس مواد تشکیل‌دهنده رنبه‌بندی شدند. 


حدول ۲- واحدهای سنگ‌شناسی منطقه موردمطالعه 


ردیف واحد سنگ‌شناسی ردیف واحد سنگ‌شناسی 
مخروطافکنه‌ها و تراس‌های آبرفتی کم مارن؛ ماسه‌سنگ آهکی» Sal Keer‏ 
Mm,s,l ۱۱ ۲ Qf \‏ 7 
ارتفاع جدید کنگلومرا 
TRJs ۲‏ ماسه‌سنگ و شيل Ek \y‏ توف و شیل توفی 
مخروط افکنه‌ها و تراس‌های آبرفتی شیل Sal‏ شیل سیلتی با میان al‏ توف 
Eksh ۳ 5 Qftı y‏ ۳0 
مرتفع قدیمی شیشه‌ای 
Qcf 3‏ کفه رسی ۱ PeEz‏ سنگ‌آهک ریفی و مارن ژییسی 
Mur ۵‏ مارن؛ مارن گچ‌دار: ماسه‌سنگ و کنگلومرا ۱۵ Sal Kw Ktzl‏ اورییتولین‌دار 
Murmg 1‏ لایه‌های مارن گچ دار ۳ Pr‏ سنگ‌آهک 
Ebv ۷‏ سنگ‌های Ekgy W ALL‏ ژیپس 
دولومیت و Sal Kee‏ ضخیم لایه با میان آهک تریلوبیت‌دان دولومیت. ماسه‌سنگ و 
COm \A pC-Cs A‏ 
Ek.a q‏ شیل آهکی با میان لایه‌های توف ۱۹ pCk‏ شیل اسلیتی و ماسه‌سنگ کوارتریتی 
Eat ۱۰‏ توف آندزیتی 


کاربری اراضی: استفاده از زمین برای هدف‌های مختلف صورت می‌گیرد و تاثیر هر یک از این عوامل بر روی 
فرونشست می‌تواند متنوع باشد. کاربری‌های کشاورزی و شهری یکی از عوامل اصلی در افزایش برداشت آب‌های 
زیرزمینی و بروز فرونشست است. همچنین وزن ناشی از این سازه‌های صنعتی سبب تراکم لایه‌های خاک و نشست 
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سطح زمین می‌شود (ژو" و همکاران» ۲۰۲۰؛ ندیری و همکاران ۱۳۹۷). لذا نقشه کاربری اراضی در مقیاس 
۰ از سازمان جنگل‌هاء مراتع و آبخیزداری ایران گرفته شد. همچنین کاربری کشاورزی» مناطق مسکونی و 
مناطق صنعتی از آن استخراج گردید. 

عمق سنگ‌بستر: در منطقه‌ای که سنگ‌بستر در عمق کمی قرار دارده به دلیل ضخامت کم آبرفت. حفر چاه 
امکان‌پذیر نیست. بر این اساس» آب‌های زیرزمینی به مناطقی منتقل می‌شوند که ضخامت بیشتری دارند که منجر به 
افزایش فرونشست و جابه‌جایی افقی Leal)‏ در آن منطقه می‌شود (پژوهشکده مطالعات و تحقیقات منابع آب AYAY‏ 
در این پژوهش لاه عمق سنگبستر از نقشه جهانی عمق سنگ‌بستر 
(https://data.isric.org/geonetwork/srv/api/records/f361 17ea-9be5-4afd-bb7d-7a3e77bf392a)‏ 
تهیه گردید. مقادیر ان پیوسته و بین صفر تا ۲۳۱ متر است. 

فاصله از جاده: محورهای ارتباطی با وجود عبور و مرورهای زیاد با ایجاد امواجی به طرفین و عمق زمین, به 
فروپاشی ذرات خاک در محل‌هایی که رطوبت کاهش پیدا کرده است ایجاد فرونشست زمین و فروجاله‌ها را تقویت 
می کند (عرب عامری" و همکاران ۲۰۲۱؛ حکیم " و همکاران. ۲۰۲۰). نقشه جاده‌ها از نقشه توپوگرافی رقومی شده 
۰ منطقه موردمطالعه استخراج و سپس با استفاده از ابزار Distance‏ در نرم‌افزار ArcGIS‏ نقشه فاصله از 
جاده‌ها تهیه گردید. 

تراکم چاه‌ها: میزان تراکم چاه‌ها و قنات‌ها با برداشت بی‌رویه از آب‌های زیرزمینی و میزان فرونشست ارتباط 
مستقیم دارد (جینه " و همکاران, ۹ در این پژوهش, موقعیت مکانی چاه‌ها به‌صورت »و لاز شرکت آب 
منطقه‌ای استان البرز (۱۳۹۹) دریافت گردید و پس از ورود به نرم‌افزار ArcGIS‏ میزان تراکم چاه‌ها بر اساس تعداد 
حلقه چاه در کیلومترمربع محاسبه گردید که مقادیر آن بین صفر تا ۲۲۳ حلقه چاه در کیلومترمربع است. 

تراکم زهکشی: میزان تراکم زهکشی با میزان فرونشست ارتباط مستقیم دارد (حکیم و همکاران ۲۰۲۰). تراکم 
زهکشی يا به عبارتی نسبت طول آبراهه‌ها به مساحت حوضه با استفاده از نقشه رقومی» ارتفاع و اکستنشن 
Archydro‏ در نرم‌افزار AreGIS‏ محاسبه گردید که مقادیر آن بین ۰/۰۶ تا ۲/۱۶ کیلومتر در کیلومترمربع است. 

ضخامت آبخوان: ازآنجایی که ضخامت آبرفت فاصله بین سطح زمین تا سنگ کف است و ضخامت آبخوان 
فاصله بین سطح ایستابی تا سنگ کف. لذا با کم کردن عمق سطح ایستابی از ضخامت آبرفت. ضخامت آبخوان را 


می‌توان برآورد نمود. به عبارتی» بخشی از ضخامت آبرفت را که اشباع از آب هست. می‌توان ضخامت آبخوان در 
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نظر گرفت. در این پژوهش, به دلیل موجود بودن نقشه نسبتاً دفیق ضخامت آبرفت و عمق سطح ایستابی, این روش 
مورداستفاده قرار گرفت. در نقشه پیوسته ضخامت آبخوان, مقادیر بین صفر تا ۲۸۵ متر است. 

عمق سطح آب زیرزمینی: میان عمق سطح آب زیرزمینی و میزان فرونشست ارتباط مستقیم وجود دارد. نقشه 
نقاط هم‌عمق پیزومترها از شرکت آب منطقه‌ای استان البرز (۱۳۹۹) دریافت گردید و با استفاده از روش IDW‏ نقشه 
پیوسته عمق سطح آب زیرزمینی تهیه شد که مقادیر آن بین ۱۰ تا ۱۶۵ متر است. 

۳-۲-۲- وزن‌دهی معیارها بر اساس روش BWM‏ 


روش بهترین- بدترین نخستین بار در سال ۲۰۱۵ توسط رضایی ارائه گردید. این روش برمبنای مقایسات زوجی 
و بهره‌گیری از مدل برنامه‌ریزی خطی استوار شده است. در این روش بهترین و بدترین شاخص توسط تصمیم 
گیرنده مشخص می‌شود و مقایسه زوجی بین هر یک از این دو شاخص (بهترین و بدترین) و دیگر شاخص‌ها 
صورت می‌گیرد. همچنین در این روش یک فرمول برای محاسبه نرخ ناسازگاری جهت بررسی اعتبار مقایسه‌ها در 
نظر گرفته شده است. مراحل این روش به شرح ذیل است (رضایی  (TAT‏ 

گام ۱: مجموعه زیرمعیارهای تصمیم‌گیری تعیین می‌شود. 

گام ۲: بهترین و بدترین زیرمعیار با اعداد ۱ تا ٩‏ مشخص می‌شود. 

گام T‏ ارجحیت بهترین معیار نسبت به دیگر معیارها به‌صورت Ap = (Api. ۵82 .... Apn)‏ مشخص می‌شود. 


گام RA‏ ارححیت بدترین معیارها نسبت به دیگر معیارها بهصورت Aw = (aw1. aw2 .... Awn)!‏ مشخص 


می‌شود. 
گام ۵: تشکیل بردار وزن بهینه به‌طوری aS‏ 
Wp Wj‏ 5 
mInmax; {ee = apj . Ww = awl} 0‏ 
st: 4.‏ 
j Wad:‏ 


W; = 0. for all j 

همچنین. مسئله را می‌توان به‌صورت زیر بازنویسی کرد. 
min č‏ 
s:t:‏ 
WB‏ 


Wj 


< €. for all j (Y) 


Ww 


W; 
h- ew > § forall j < ۷ = 1; W; > 0. forall j 
7 
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۲-۲-- اعتبارسنحی مدل 

به‌منظور اعتبارسنجی از منحنی تشخیص عملکرد نسبی (ROC)‏ مساحت سطح زیر این منحنی (AUC)‏ استفاده 
شد. محور عمودی منحنی نشان‌دهنده درصد L Sensitivity‏ حساسیت است و محور ol Ël‏ نشان‌دهنده 
Specificity‏ يا ویژگی. مقدار AUC‏ نشان‌دهنده احتمال این موضوع است که یک پیکسل انتخاب شده به‌طور 


دقت از نقشه‌های بهنه‌بندی agi‏ شده است (پرادهان" و همکاران, ۲۰۱۶). 


۳- بحث و نتایج 


در این پژوهش. ۲ بعد فیزیکی- شیمیایی و اجتماعی- اقتصادی مشتمل بر ۵ معیار و ۱۹ زیرمعیار برای تعیین 
پهنه‌های مستعد فرونشست تدوین گردید و در ادامه بر اساس روش Z-SCOTE‏ استاندارد شدند. شکل‌های (۳) تا 
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شکل -i‏ طبقهب: 5 
شکل ۶- طبقه‌بندی معیارهای هوا و اقلیم در پهنه‌بندی فرونشست 
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40000 
3 0 


شکل 


طبقه‌بندی معیارهای کاربری اراضی در پهنه‌بندی فرونشست 


شکل ۷- طبقه 


بندی معیارهای زمین‌شناسی در پهنه‌بندی فرونشست 


در ادامه ۲۰ عدد پرسشنامه بین کارشناسان ژئومورفولوژی و زمین‌شناسی توزیع و با استفاده از مدل BWM‏ 
وزن معیارها و زیرمعیارها محاسبه گردید. برای این کار در ابتدا بر اساس نظر کارشناسان مهم‌ترین و کم اهمیت‌ترین 
زیرمعیار مشخص شد؛ در ادامه پرسشنامه‌ها در دو بخش (بخش اول مقایسه بهترین زیرمعیار با سایر زیرمعیارها و 
بخش دوم مقایسه ple‏ زیرمعیارها با بدترین زیرمعیار) بین کارشناسان توزیع شدند تا بر اساس اصطلاحات زبانی آن 
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شد تا بر اساس OF‏ وزن زیرمعیارها محاسبه گردد. سپس با ميانگین گیری از تمام خبرگان» ضریب اهمیت و یا وزن 
هر یک از معیارها Gilles‏ جدول (۳) به‌دست آمد. از آنجایی که "2 بیانگر میزان سازگاری یک مقایسه می‌باشد. این 
شاخص برای همه کارشناسان در این پژوهش نزدیک به صفر است که می‌توان نتیجه گرفت نتایج به‌دست آمده 


سا زگاری بسیار خوب داشته است و قابل قبول بوده است. 


جدول ۳- میزان اهمیت معیارهای مؤثر بر فرونشست دشت هشتگرد 


بعد معیار زیرمعیار وزن وزن نهایی رتبه 

افت سطح آب زیرزمینی AYI‏ ۷۹ ۳ 

برداشت آب زیرزمینی ۳۷۵ Ak‏ ۱ 

ضخامت آبخوان AY‏ ۰/1۸ ۸ 

هیدرولوژی ضخامت آبرفت ۲ y N‏ 

\s Yt 0۸ تغذیه آبخوان‎ YM 

4 t ۰/۰۹۸ تراکم چاه‌ها‎ oe ۲ 

۳ AL iY نی تراکم زهکشی‎ ds 
۳ ۰/۳۵5۸ 

عمق اب زیرزمینی YAVA AY‏ 3 

زمین‌شناسی نوع سازند زمین‌شناسی YNY‏ ۰/۱۱0۷ ۲ 

۱۱ YY W/W عمق سنگ‌بستر‎ ۲ 

میزان بارش سالانه EY‏ 4 4 

متوسط دمای سالانه ۰/۱ ۹ ۱۷ 

\o AY ۱0۷ نوع اقلیم سرزمین‎ l 

\y JAN “yyy کاربری کشاورزی‎ 

۳ ie At کاربری شهری‎ eee 

اجتماعی- اقتصادی i‏ کاربری صنعتی WV WT‏ ۵ 

۱۸ ۹ ۵ فاصله از رودخانه‌ها‎ EY 

۱۹ Yeg /\YO فاصله از جاده‌ها‎ pr 

۱ ۰۸ ۰/۵ فاصله از گسل‌ها‎ a 


بر اساس روش BWM‏ و با توجه به نظر کارشناسان» عوامل برداشت از آب زیرزمینی CYNA)‏ نوع سازند 
زمین‌شناسی (۰/۱۵۷) افت سطح آب زیرزمینی (۰/0۷۹) و عمق آب زیرزمینی (VA)‏ اهمیت بیشتری بر پتانسیل 
خطر وقوع فرونشست در دشت هشتگرد دارند. در این مرحله» لایه‌ها با استفاده از شاخص همپوشانی وزین جهعت 
ag‏ نقشه پهنه‌بندی فرونشست ادغام شدند که در شکل (A)‏ نشان داده شده است. 


سال یازدهم پهنه بندی خطر وقوع فرونشست زمین در دشت ... ۱:۳ 


460000 480000 500000 


خبلی کم ۷ 
کہ" 
مثوسط ° 
زباد * 


خیلی زیاد * ۳۸/۵۱ 


با سس 
13 32565 0 
Kilometers‏ 


500000 480000 460000 440000 
شکل ۸- پهنه‌بندی پتانسیل فرونشست بر اساس روش شاخص هم‌پوشانی وزنی در دشت هشتگرد 


با توجه به بررسی‌های صورت MVS‏ درصد از مساحت منطقه موردمطالعه به طبقه خیلی کم‌خطر 
۱ درصد به aib‏ خطر کم ۳۱/٤۹‏ درصد به aib‏ خطر متوسط. ۱۱/۱۱ درصد به طبقه پرخطر و در نهایت ۷/۱۹ 
درصد به طبقه خیلی پرخطر اختصاص بافته است. بر Gb‏ نقشه خروجی. نیمه مرکزی دشت هشتگرد حساسیت 
بیشتری برای فرونشست دارد و امکان وقوع این مخاطره در این بخش حوضه بسیار زیاد می‌باشد. هم‌پوشانی این 
نقشه با نقشه معیارها نشان داد که بیشتر مناطق آسیب‌پذیر به فرونشست در نزدیکی گسل‌هاء مناطقی با بیشترین 
برداشت آب زیرزمینی» اراضی کشاورزی و در مخروط افکنه‌ها و تراس‌های آبرفتی مرتفع قدیمی اتفاق افتاده است. 
درنهایت به‌منظور ارزیابی دقت و کارایی مدل برای پهنه‌بندی خطر فرونشست. از منحنی ROC‏ و سطح زیر منحنی 
AUC‏ استفاده شد. بدین منظور نقاط برداشت شده و cote polie‏ متناظر با آن‌ها در نقشه نهابی پتانسیل فرونشست 
وارد محیط نرم‌افزار SPSS‏ شد و منحنی ROC‏ محاسبه شد. سطح زیر منحنی برابر با ۰/۹ درصد است که عملکرد 
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AUC of‏ سس 
WOI=0.9‏ 


Reference line 


5 7 9 111315171921 23252731333537 4143454749 51 
False Positive 6 


شکل -٩‏ نمودار ارزیابی عملکرد مدل بر اساس منحنی‌های ROC‏ 


6- نتیجه گیری 

در تحقیق حاضر برای ارزیابی پتانسیل وقوع فرونشست در دشت هشتگرد از ۱٩‏ زیرمعیار تأثیرگذار در وقوع 
فرونشست. شامل میزان افت آب زیرزمینی» برداشت آب زیرزمینی» ضخامت آبخوان. ضخامت آبرفت تغذیه 
آبخوان» تراکم daela‏ تراکم زهکشی» عمق سطح آب زیرزمینی» نوع سازند زمین‌شناسیء عمق سنگبستر میزان 
بارش سالانه» متوسط دمای سالانه نوع اقلیم سرزمین» کاربری کشاورزی» کاربری شهری» کاربری صنعتی» فاصله از 
آبراهه» فاصله از جاده و فاصله از گسل استفاده و با مدل BWM‏ وزن‌دهی شدند. بر اساس مدل BWM‏ و با توجه به 
نظر کارشناسان» عوامل برداشت از آب زیرزمینی (۰/۲۱۹ نوع سازند زمین‌شناسی (۰/۱۵۷» افت سطح آب 
زیرزمینی (CVA)‏ و عمق آب زیرزمینی OVA)‏ اهمیت بیشتری بر پتانسیل رخداد فرونشست در دشت هشتگرد 
دارند. 

درنهایت» با تلفیق وزن معیارها و نقشه‌های استاندارد شده معیارهاء نقشه نهایی پتانسیل مناطق مستعد وقوع 
فرونشست در دشت هشتگرد به‌دست آمد. بر اساس نتایج به‌دست آمده مناطق مرکزی دشت هشتگرد ریسک بالایی 
aaa‏ اک dracon‏ مین پوس مات هر یک او گنها als OL‏ ۱/۷۲ دراد 
مساحت منطقه موردمطالعه به طبقه خیلی کم خطر ۳۸/۵۱ درصد به طبقه خطر کې ۳۱/4٩‏ درصد به طبقه حطر 
متوسط. ۱۱/۷۱ درصد به طبقه پرحطر و درنهایت ۷/۱۹ درصد به طبقه خیلی پر خطر اختصاص adh‏ است. بنابراین» 
بر اساس نتایج نقشه حساسیت تهیه شده فرونشست‌ها و نیز با توجه به تعیین میزان اهمیت عوامل موثر در وقوع 
فرونشست. نتایج این پژوهش می‌تواند به برنامه‌ریزان و تصمیم‌گیران در پیشگیری از این مخاطره محیطی و 
برنامه‌ریزی برای مدیریت منابع آب و جلوگیری از برداشت بی‌رویه از آب‌های زیرزمینی کمک شایانی نماید. 


سال یازدهم پهنه بندی خطر وقوع فرونشست زمین در دشت aus‏ ۱1۵ 


همچنین, نگارندگان بر این باورند که در تحقیقات آینده» استفاده از روش‌های دیگر همچون ماشین یادگیری و 
پایش و بررسی جامع وضعیت فرونشست منطقه با بهره‌گیری از AUT‏ سری‌های زمانی تصاویر راداری می‌تواند کمک 


شایانی در طرح‌های آمایش سرزمین و طرح‌های آبخیزداری داشته باشد. 
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